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Abstract 

We study the dynamics of the "externally" forced and damped Fermi-Pasta-Ulam 
(FPU) ID lattice. The forcing has the spatial symmetry of the Fourier mode with 
wavenumber p and oscillates sinusoidally in time with the frequency tu. When to is 
in the phonon band, the p-mode becomes modulationally unstable above a critical 
forcing, which we determine analytically in terms of the parameters of the system. 
For uj above the phonon band, the instability of the p-mode leads to the formation 
of a travelling multibreather, that, in the low-amplitude limit could be described in 
terms of soliton solutions of a suitable driven-damped nonlinear Schrodinger (NLS) 
equation. Similar mechanisms of instability could show up in easy-axis magnetic 
structures, that are governed by such NLS equations. 
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1 Introduction 



Chains of anharmonic oscillators constitute an ideal testing ground for the 
different theoretical approaches to the study of the dynamics of many-degrees- 
of-freedom systems. For instance, the first numerical study of the approach to 
equilibrium was performed by Fermi, Pasta and Ulam (FPU), initializing the 
chain in a long wave-length state and looking at mode energy sharing at long- 
times [1] . Later on, the role of certain exact nonlinear standing and travelling 
wave solutions of such Hamiltonian lattices was discussed in relation with 
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chain excitations at shorter wavelengths [2-4]. This led to the discovery of 
chaotic breathers [5], spatially highly localized and erratically moving objects, 
which appear in the process of thermalization at small energy densities (these 
were in fact previously found in Ref. [6] without emphasizing too much their 
role for the approach to equilibrium). This discovery motivated a series of pa- 
pers [7-9] where the formation, lifetimes and destruction of these "intrinsically 
localized" structures were studied in full detail. Very recently, a different ap- 
proach was followed in which, instead of looking at the fate of a specific initial 
state of such FPU lattices, e.g. the highest frequency zone-boundary Fourier 
7r-mode, a driving term with this spatial symmetry was directly added to the 
Hamiltonian equations of motion [10]. Damping was included to allow the for- 
mation of stationary states and it was taken to be weak to remain close to the 
Hamiltonian case. As a result the zone-boundary mode was found to be stable 
and locked to the driving field below a certain critical forcing, which is indeed 
analytically calculable. Above such forcing, a standing modulated wave forms 
for driving frequencies that are below the band edge, while a static multi- 
breather (a train of breathers) develops for forcing frequencies above the band 
edge. The patterns emerging from the instability of the zone-boundary mode 
are always almost static, because the group velocity of the carrier wave is zero. 
They can acquire a nonzero velocity due to other nonlinear mechanisms [7] . 

In this paper we study a damped FPU ID lattice forced by a travelling wave 
field with wavenumber p ^ n. An instability mechanism similar to the one 
mentioned above leads to the formation of travelling modulated waves for in- 
band frequencies. When the forcing frequency is instead above the phonon 
band, a moving multibreather arises after instability of the p-mode. This state 
can indeed be thought as a moving envelope soliton train (the carrier wave 
having the symmetry of the p-mode) , since the envelope of these "intrinsically 
localized" moving objects can be put into a close relation with the soliton 
solutions of a suitable externally driven and damped nonlinear Schrodinger 
(NLS) equation. Such type of equation arises in the small amplitude limit of 
the ac driven damped sine-Gordon (SG) system; it is discussed in the literature 
and its analytical and numerical solutions are known [11-13]. 

In Section 2 we introduce a Fourier mode representation of the FPU lattice 
and we study it in the rotating wave approximation and for weak damping, 
obtaining an explicit analytical expression for the critical amplitude f cr in 
terms of the parameters of the system. Section 3 is devoted to the description 
of the travelling multibreather, both numerically and, analytically, in terms of 
a suitable NLS equation. In Section 4 we draw some conclusions. 
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2 Stability of the travelling wave 



The equations of motion of the "externally driven" damped FPU chain read 
as follows: 



ii n = U n+ i + W n _i - 2u n + (u n+1 - U n f +(u„_i - U n f 

-iu n + fcos(ut + pn), (1) 

where u n is the displacement of n-th oscillator with respect to its equilibrium 
position. Periodic boundary conditions, u n+ N = u n , are assumed, with N being 
the number of oscillators. Dimensionless units are used such that the masses, 
the linear and nonlinear force constants and the lattice spacing are taken equal 
to unity. The forcing and damping strengths are gauged by the parameters / 
and 7, respectively; u and p are the driving frequency and wavenumber. For 
technical reasons, to be discussed below, we restrict to tt/2 < \p\ < n. As in 
Ref. [10] we use the following expansion: 



Z k 



\ E [a* e i{Mn) + at k e"*^"*")] , (2) 



where are complex mode amplitudes and — 7r < k < tt the corresponding 
wavenumbers. The equations of motion for the amplitudes a^'s are obtained by 
substituting Eq. (2) into Eq. (1). Similarly to what is done for the undamped 
case [14], a considerable simplification is achieved by neglecting higher-order 
harmonics that are produced by the cubic force terms. This is the so-called 
rotating wave approximation, that is valid for small amplitudes, \a k \ 1. 
Moreover, in the limit of weak damping 7 < w and neglecting inertial terms 
dk, we find the following set of approximate equations: 

- 2iud k - iu^a k = (u 2 k - u 2 )a k + f5 k , p + 6 £ G k quQ2 a qi a q2 a+ +q2 _ k , (3) 

where 

u 2 k = 2(1 -cos k), (4) 
5 k>p is equal to one for k = p and zero otherwise and 

^91,92 = ^t 1 + cos (^i + Sa) + cos(£; - q 2 ) + cos(A; - q x ) 

— cos k — cos qi — cos q 2 — cos(k — q 1 — q 2 )\ . (5) 
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Let us begin by considering solutions where only the p-mode is excited. Such 
solutions, corresponding to nonlinear travelling waves with velocity v = oo/p, 
are numerically observed to exist and to be stable below a certain critical 
forcing f cr . Indeed, when we initialize the lattice at equilibrium and we impose 
random velocities, all modes with k ^ p damp out on a time-scale set by the 
value of 7 _1 while \a p \ rapidly grows and finally approaches a constant value. In 
this asymptotic regime the system is well described by the "internal dynamics" 
of the p-mode, given by the following equation: 

3 

— 2iuoa p — ioj^dp = (ujp — uo 2 )a p + f + -oo p a p \a p \ 2 . (6) 



The asymptotic amplitude a p is then obtained by solving for the stationary 
solution of Eq. (6), which originates the algebraic equation below: 



which can be numerically solved for both the modulus and the phase of a p . 
In particular the squared modulus of the solution of the cubic 

equation: 



9 o o 3 



-u~z~ - y v {u? - uj 2 p )z 2 + [(u 2 - ^ p f + 7 z = p . (8) 



In reality, also the higher odd harmonics of p will be excited, but they can be 
neglected when \a p \ 1. One can easily ascertain, as it also happens for the 
p = it case [10], that Eq. (8) admits a single real root only for |u;| < cj* where 



a;* ~ lu. 



' P + >/3t/2, (9) 



while three distinct roots may otherwise exist. We will not analyse in this 
short note what happens when three solutions of the consistency equation (7) 
coexist for a given parameter set. It turns out that, the destabilization of the 
p-mode can be described by considering the lowest amplitude solution of (6). 

The p-mode nonlinear solution is stable for / < f cr (uj^,p). It is possible to 
get an analytic expression for / cr -(^,7,p) in the range < u < u*. This is 
accomplished by solving the following set of equations, that one obtains when 
linearizing Eqs. (3) around the stationary p-mode solutions (i.e. neglecting all 
mode interaction terms which are nonlinear in the ajt's) 



-2iuja k - iu>ya k = (u 2 k - u 2 )a k + -u 2 p B Kp a 2 p al p _ k 
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where 



3 

-2iuja+ p _ k - i^a% p _ k = -{u%_ k - uj 2 )a+ p _ k - -u 2 p B k:P (a+) 2 a k , (10) 



3 

^1 = (! + 2^pl a p! 2 )^' B k,v = cos ( k -P)~ cos P ( n ) 



are the frequency of the fc-th mode shifted by the interaction with the j9-mode 
and B ktP derive from G qiq2 in Eq. (5), when the interaction is restricted to 
the mode subset k,p, 2p — k. Since B k , p = B 2p -k, P , the two Eqs. (10) can be 
obtained one from the other for k — > (2p — k), hence they can be solved looking 
for symmetric solutions of the form a k ~ a 2p -k ~ exp(u k t). The relevant 
branch of the eigenvalue spectrum reads 



7 , - w l ~ vip-k 
Vk = ~2 +l Au 



+ 2M V 4^ |Gp|4 " ~ 2°** + ^- k) ) ■ (12) 

The growth rate Re{u k } is maximal when the square root in the above ex- 
pression attains its maximum value, i.e. when the "resonance" condition 

2u 2 = u 2 k + u 2 2p _ k (13) 



holds. The latter, together with the expression for (b\ in (11), allows to get 
the following equation for the wavenumber k* of the most unstable mode 



u 2 



COS (p-K) = I 1 O _|_ Q 21 12 J > ( M ) 

cosp \ 2 + 3iu^\a p \ 2 ) J 



in terms of the amplitude of the p-mode. Moreover, if we let i?e{^} = in 
Eq. (12) and use the definition of B kp in (11), we obtain the following equation 
for the squared amplitude of the p-mode: 



3u p \ cos(&;* — p) — cosp 



Solving the set of Eqs. (14) and (15) for both k* and \a p \ one obtains the 
critical amplitude \a p \ cr of the p-mode, above which modulational instability 
occurs. Finally, replacing the expression for z cr = \a p \ 2 r in Eq. (8), we obtain 
fcr{li ljJ iP)i as announced above. We do not display its explicit expression 
because of the rather bulky form. All this is of course valid when only one 
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solution of the cubic equation (8) is present, i.e. for uj < cj*. Moreover, Eq. (14) 
has a solution in this full u range only when tt/2 < \p\ < it, as can be easily 
checked. In this range of parameters we can get rather general expressions both 
for the stable p-mode pattern and for the expression of the critical forcing f cr . 

When / > f cr the p-mode is modulationally unstable and, presumably, as 
for the p = 7i case [10], the system saturates asymptotically, due to nonlinear 
effects, into a state where the triplet of modes k*,p, 2p — k* is excited, resulting 
into a more complex travelling modulated wave. Its expression can possibly be 
computed, solving for the stationary solution of the coupled equations for this 
triplet of modes. This calculation is out of the scope of the present short note. 



3 Travelling mult ibr eat hers 



As it was mentioned above, in the range \uj\ > cu* = iv p + y/3^/2, three roots of 
Eq. (8) are present and the instability process is determined by the "internal" 
dynamics of the p-mode, as defined by Eq. (6). After instability, several modes 
grow at the same time and no periodic pattern develops. However, numerical 
simulations show that a travelling multibreather can form just above threshold. 
An example is shown in Fig. 1, where we plot the local energy 



1 -2 1 

K = - 2 u n + - 



u n+l ~~ u r 



+ 



K+l - U n ) 4 + (U n - Un-xf 



(16) 



vs. the lattice position sampled at the period of the forcing. The corresponding 
spatial Fourier spectrum is shown in Fig. 2. The broad band structure of the 
spectrum reflect the non perfect periodic arrangement of the localized peaks 
in Fig. 1. 

Such states can be described in terms of soliton solutions of an associated 
suitable driven-damped nonlinear Schrodinger (NLS) equation. Let us first 
make the following definition: 



1 



u, 



v (n,t)e l ^ t+pn) +a+(n,t)e 



-i(tjjt+pn) 



(17) 



where a p and its conjugate are smooth functions of n. Such an assumption is 
possible if the wavepacket is concentrated around the driving mode Ak <^ ir. 
Substituting (17) into the equations of motion (1) one gets: 
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w 2 <9 2 a p 3 4 , 



After performing the following re-scalings: 





^ 2 -wj 




[3 






'uo 2 — 



uj 2 — UJ 2 



n — vt), 



a; 



c<-> 2 — a; 2 . 



7, 



8 (w 2 - w 2 ) 3 / 2 



and choosing a reference frame moving with velocity v = du p /dp = sinp/uj p , 
Eq.(18) reduces to the well studied "externally" driven (or ac driven) damped 
NLS equation [11,12]: 



(19) 



Exact soliton solutions of this equation can be obtained for 7' = 0, see Eqs.(37- 
40) of Ref. [11]. Moreover, multisoliton solutions are also derived in Ref. [12]. 
What we observe in Fig. 1 might well be a superposition of such solutions to 
form a train of "intrinsically localized" structures. However, one should bear 
in mind that NLS solutions can describe only low amplitude states. Therefore, 
they can be only a fair approximation of the pattern displayed in Fig. 1, which 
shows high amplitude localized peaks. 

In Fig. 3 we plot the speed of the travelling multibreather as a function of 
the wavenumber of the forcing p-mode, which compares well with the group 
velocity of the corresponding linear waves, showing that nonlinear effects are 
negligible in this parameter range. 



4 Conclusions 



We have shown how to extend the analysis of the 7r-mode solution in an "ex- 
ternally" forced and damped FPU ID lattice performed in Ref. [10] to that 
of a generic p-mode solution with 7r/2 < \p\ < ir. Since the carrier wave has a 
nonzero velocity, we find, for weak forcing, solutions corresponding to nonlin- 
ear travelling waves. These solutions become unstable above a critical forcing 
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f cr , of which an analytical expression can be derived in terms of the param- 
eters of the system: the frequency of the forcing uj, the damping constant 7 
and the wavenumber p. Above f cr the p-mode solution is modulationally un- 
stable and the system generates complex patterns after nonlinear saturation: 
travelling modulated waves for in-band frequency forcing uj < cj* or travelling 
multibreathers for out-band forcing u > uj*. We suggest that the multibreather 
pattern could be described by the soliton solutions of a driven-damped non- 
linear Schrddinger equation. Travelling soliton trains have been also recently 
observed for the weakly damped parametrically driven NLS equation [15]. Our 
study also hints at the existence of similar mechanisms for the generation of 
moving breathers in more complex physical systems, which are known to be 
governed by the same driven-damped NLS or Sine-Gordon equation, e.g. one 
dimensional easy-axis magnetic structures [16]. 
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Driven FPU f=0.075 omega=2.05 gamma=0.1 k=2.4542 




lattice site 



Fig. 1. Travelling multibreather pattern generated after the modulational instability 
of the p-mode for a lattice of N = 512 sites. Here, u = 2.05, / = 0.075, p = 2.4542. 

4.0 i— — — — ^ — — 1 



2.0 




Fig. 2. Spatial spectrum of the travelling multibreather. = \Uk\ + ^ 2 |^fe|) where 
Uk is the fc-th component of the Fourier spectrum of the displacement field u n . Same 
parameters as in Fig. 1. 
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Fig. 3. Velocity of travelling multibreathers versus the wavenumber of the forcing. 
The solid line is the group velocity of linear waves. 
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